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Abstract 

We have investigated the photoneutron cross section of the isotopes 148 > 150 Nd, 154 Sm, and 
,160 Gd close to the neutron emission threshold in photoactivation experiments at the Darm- 
stadt superconducting electron linear accelerator S-DALINAC. Naturally composed targets were 
activated with a high-intensity bremsstrahlung beam at various energies and the reaction yields 
have been determined by measuring the activity of the produced radioactive isotopes with HPGe 
detectors. The results are compared to two different statistical model calculations. 
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I. INTRODUCTION 



The nucleosynthesis of the elements heavier than iron is dominated by the slow neutron 
capture process (s process) and the rapid neutron capture process (r process). However, 
35 proton-rich stable isotopes cannot be produced in either of these two processes and 
are believed to be mainly synthesized by a combination of photodisintegration reactions, 
i. e. (7, n), (7,p) and (7,0) reactions, in the explosive scenario of the p process |2J. An 
accurate description of the nucleosynthesis within these processes demands reliable input 
from nuclear physics. 

Calculations of astrophysical reaction networks within the p process have to account for 
a huge number of reaction rates. Many of the isotopes involved in the network are off 
the valley of stability and, therefore, are not accessible by experiments in many cases. 
But even for the stable nuclei experimental data for photodisintegration reactions in the 
astrophysically relevant energy region, i. e. close to the particle emission threshold, are 
rare. Hence, the reaction rates mainly need to be adopted from theoretical predictions. 
It is mandatory to test the predictive power of the calculations for those isotopes, which 
can be studied in the laboratory, and to prove the reliability of the predictions for the 
extrapolation to unstable nuclei. 

The situation is much better for the s process. Since for the unstable isotopes within the 
s-process reaction network the /3-decay rate usually highly exceeds the neutron capture 
rate, the reaction path follows the valley of stability and, hence, mainly stable isotopes 
are involved. Therefore, extensive experimental studies have been carried out in neutron 



'or a wide range 
J]. With this 



capture experiments in the last decades and neutron capture cross sections 
of isotopes have become available with uncertainties of only a few percent 
enormous amount of nuclear physics input the understanding of s-process nucleosynthesis 
and its astrophysical sites has largely improved. However, for some unstable isotopes 
along the reaction path, the so-called branching points, /3-decay and neutron capture 
become competetive. Since the branching ratio of /3-decay and neutron capture is highly 
sensitive to the physical conditions, e. g. temperature and neutron density, in the stellar 
environment during the s process, these branching points serve as an excellent test for 
s-process nucleosynthesis models. Thus, it is essential to determine the neutron capture 
cross sections of these isotopes with high accuracy. Unfortunately, due to their short half-life 
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neutron capture experiments can hardly be performed and, hence, further improvements 
of nucleosynthesis models are still hampered by the relatively large uncertainties of the 
cross section predictions for these isotopes 5|. Nevertheless, information for the neutron 
capture cross sections of branching points can be derived from studying the (7, n) reaction 
of the stable neighbouring nucleus. The idea to determine (n, 7) reaction rates via 
photodisintegration experiments has already been presented in jf]]. 

Experiments with real photons provide a well-suited tool to study photodisintegration cross 
sections and ground-state reaction rates. In previous experiments we have concentrated our 
investigations on isotopes of mass A > 186 and confirmed the validity of several theoretical 
predictions in this mass region , llO| . Average deviations between experiment and 



theory were typically less than 30%. 

In this paper we want to address the photodisintegration reactions in the rare earth region 
(see Fig. [I]). This region is of particular interest for the s process, because the isotopic 
abundance pattern is shaped by several branchings at the unstable isotopes of Nd, Pm, 
Sm, Eu and Gd. Since the relative abundances of the involved isotopes are known with an 
uncertainty of better than two percent HI], this mass region is exceptionally suited to test 
the stellar nucleosynthesis models with high accuracy. To provide precise nuclear physics 
input, neutron capture cross sections of several branching points have been measured with 



high accuracy in the last years 
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But still one has to rely on theoretical 



predictions of the cross sections of short-lived branching points, e. g. 147 Nd and 153 Gd. 
Therefore, we have investigated the inverse photoneutron reactions of 148 Nd and 154 Gd in 
order to improve the reliability of these predictions. 

Furthermore, the abundances of some isotopes in this region receive large contributions 
from the p process, e. g. a 33% contribution is predicted for 152 Gd ll|. For a comprehensive 
description of the nucleosynthesis in this mass region these contributions cannot be 
neglected and need to be studied in detail. Hence, the aim of this work is also to provide 
reliable data for the photoneutron reaction rates of some selected isotopes to prove the 
predictive power of the theoretical calculations in this mass region. 



We have performed photoactivation experiments to study the (7, n) reactions of the iso- 
topes 148 ' 150 Nd, 154 Sm and 154,160 Gd in the astrophysically relevant energy region close to the 
neutron emission threshold at the superconducting Darmstadt electron linear accelerator 
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FIG. 1: The s-process flow (solid lines) in the rare earth region. The thickness of the solid lines 
indicates the strength of the reaction flow. The reaction path is influenced by several branch- 
ing points (dashed boxes). Unstable isotopes are shaded light grey (J3~) and dark grey ((3 + ,e). 
Most isotopes also have nucleosynthesis contributions from the r and p process (dotted lines), 
respectively. The photoneutron reactions of the hatched isotopes have been studied in this work. 
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S-DALINAC [16!. Our experimental setup is discussed in detail in Sec. [TTJ We explain the 
data analysis and the approach to derive astrophysically relevant ground-state reaction rates 
directly from the experiment by approximating a thermal Planck spectrum by a superposi- 
tion of different bremsstrahlung spectra in Sec. IHIl The results are presented in Sec. HVI and 
compared to two theoretical calculations based on the framework of the Hauser-Feshbach 
theory. 



II. EXPERIMENTAL SETUP 



We have irradiated naturally composed neodymium, gadolinium and samarium targets 
at the superconducting electron linear accelerator S-DALINAC. The experimental setup is 
illustrated in Fig. [2j A monoenergetic electron beam of energy E is stopped completely 
in a thick copper radiator and produces a continuous spectrum of bremsstrahlung photons 
with a maximum energy of E m3jX = Eq. The photons irradiate the target, which is mounted 
directly behind the radiator. Thin metallic discs (m = 40 — 50 mg) were used as neodymium 
and gadolinium targets, whereas samarium was available as pressed pills of SHI2O3 powder 
(m = 500 — 1750 mg). Each target had a diameter of 20 mm. The target specifications are 
summarized in Tab. U 

The photon flux intensity was determined by observing the photon scattering reaction 
11 B(7,7 / ) behind a copper collimator system with actively-shielded high-purity germanium 
(HPGe) detectors. To compare the photon flux intensity at the target position in front of 
the collimator and at the position of the 11 B(7, 7') target behind the collimator thin metallic 
rhenium discs have been irradiated simultaneously at both positions to normalize the pho ton 
flux intensity via the photodisintegration reaction 187 Re(7, n), which was studied in [10|]. A 
detailed discussion of the photon flux calibration is given in the following section. 

The advantage of irradiating the targets in front of the collimator becomes clear from 
Fig. [3] At this position the much more intense photon flux guarantees a high reaction yield 
for the irradiation even for small amounts of target material and small photodisintegration 
cross sections. 

The irradiation was performed at various energies E max starting just above the neutron 
emission threshold of 148 ' 150 Nd with E maiX = 7450 keV. The energy was increased in steps of 
150 and 200 keV up to E max = 9800 keV. The duration of each activation run was between 
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Isotope 



Form 



Abundance [%] Weight [mg] S n [keV] Analyzed 7-transition [kcV] 



/ 7 [%] 



s Nd 



metal foil 



J Nd metal foil 



154 Gd metal foil 



Sm oxide 



5.76 



5.64 



2.18 



J Gd metal foil 21.86 



22.70 



7 Re metal foil 62.60 



40-50 



7332.9 



40-50 7379.9 



40-50 8894..' 



40-50 7451.4 



500-1750 7967.6 



320-340 7363.0 



91.1 
531.0 

114.3 
155.9 
211.3 
267.7 
270.1 

97.4 
103.2 

58.0 
226.0 
348.3 
363.6 

69.7 
97.4 
103.2 

122.6 
137.2 



27.9 ± 1.1 

13.1 ±0.9 

19.2 ± 1.5 
5.93 ±0.31 
25.9 ± 1.4 
6.03 ±0.28 
10.6 ±0.5 

29.0 ±0.8 

21.1 ±0.6 

2.49 ±0.07 
0.217 ±0.002 
0.239 ±0.003 

11.8 ±0.1 

4.73 ±0.04 
0.772 ±0.019 

29.3 ±0.2 

0.603 ± 0.003 
9.47 ±0.30 



TABLE I: Specifications of targets and ca' 
intensities per decay J 7 were taken from 



ibration targets used for the activation experiments. The 
170- 



6 and 24 hours. 



A. Calibration of photon flux 

The spectral distribution of the photon flux is taken from a simulation using the Monte 



Carlo code Geant4 [18j. The absolute normalization of the photon flux intensity can be 
derived from the reaction yields of the photodisintegration reaction 187 Re(7,n) and the 
photon scattering reaction 11 B(7,7 / ), respectively. 

The reaction yield of the photoactivation using bremsstrahlung with a maximum 

energy £ max is given by: 
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FIG. 2: Schematic layout of the photoactivation setup. A monoenergetic electron beam produces 
bremsstrahlung in a massive copper radiator with a thickness of 1.2 cm (1). The photons activate 
the target of interest (x± ~ 5 cm), which is sandwiched between several calibration targets (2). 
Behind a thick copper collimator system (3) a boron target (4) is mounted (x2 ~ 150 cm) to 
monitor the photon flux via the photon scattering reaction 11 B(7,7 / ) with actively-shielded HPGe 
detectors (5). 



= N T J Nj(E, £ max ) * {j , n) (E) dE , (1) 

where Nt denotes the number of target nuclei, N^(E,E mauX ) the time- integrated photon 
flux J Q * act n^(E, E max ,t) dt for the duration of activation t act and o"( 7 ra ) is the photoneutron 
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FIG. 3: Available photon flux n 7 at our photoactivation setup for a photon energy of = 
S' n ( 148 Nd) = 7332.9 keV as a function of the maximum photon energy of the spectral distribution 
i? max . A mean electron current of I e = 25 //A on the radiator has been assumed. The photon 
flux was obtained from a simulation using the Monte Carlo code Geant4 18j and confirmed 



experimentally by the two approaches discussed in the text. The photon flux in front of the 
collimator is about 270 times more intense than behind the collimator and the energy dependence 
is nearly equal. 

cross section. Using a simulation N ljS i m (E, E mauX ) for the spectral distribution of the photon 
flux the normalization factor A^fs 1 ax for the simulation can then be derived from: 

"y ax 

jy£ max _ (7£0 (2) 

" N T j N %sim (E, £ max ) a (7 , n) (£) dE' V J 

By irradiating 187 Re both in front and behind the collimator one can determine the 
normalization for the photon flux intensity at both target positions from the well-known 
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photoneutron cross section of 187 Re and the measured reaction yield Y^ 1 ^. 

A second approach for the normalization is to observe the photon scattering reaction 
11 B(7,7 / ) with HPGe detectors behind the collimator. The reaction yield for a certain 
transition from a state of energy Ei into a state of energy Ej is given by: 



— N Bor ■ N y (Ei, E m£kX ) ■ , (3) 

where iVB or is the number of n B nuclei in the target and I^j denotes the integrated 
cross section of the observed 7-transition. With the knowledge of one can directly 
determine the photon flux intensity at various energies by observing different transitions. 
These data points can then be used for the normalization of the full photon spectrum. Figure 
H] shows the normalized photon flux distribution calculated with Geant4 in comparison to 
the experimental data points of the 11 B(7,7 / ) reaction. A mean deviation of 10% to 20% 
between simulation and experimental data was found depending on the maximum photon 
energy E max . In comparison with former simulations using Geant3 19] the shape of the 
photon flux is correctly described by Geant4 and, therefore, a correction procedure of the 
shape of the simulated photon flux distribution close to E max [ij does not need to be applied 
anymore. This proves the reliability of the simulation. 

Nevertheless, this approach only yields the normalization of the photon flux intensity 
at the target position behind the collimator, but we assumed that the same normalization 
factor was valid at the target position in front of the collimator. In addition, to confirm the 
normalization based on the 11 B(7,7 / ) reaction the photon flux intensity was derived from 
the photodisintegration reaction 187 Re(7,n) at both target positions. As seen from Fig. [5] 
both approaches are in excellent agreement. 



B. Determination of reaction yield 

In order to determine the reaction yield Y^^, i. e. the number of (7, n)-reactions occur- 
ring during the irradiation, the 7-transitions following the /3-decays of the produced unstable 
isotopes were measured offline with HPGe detectors after the activation. The detectors were 
covered by thick lead shielding to reduce natural background. The reaction yield is directly 
proportional to the peak areas of the corresponding 7-transitions in the spectra. The factor 
of proportionality can be determined with regard to the activation time, the detection effi- 
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FIG. 4: The simulated photon flux distribution at the target position behind the collimator for 
£ max = 9050 keV. The simulation was normalized to the experimental data obtained from the 
reaction yields of the photon scattering reaction 11 B(7,7') (see text for details). 

ciency and the 7-intensities. This has been discussed in detail in 

Two different setups have been used to measure the activity of the produced unstable iso- 
topes. One setup consisted of two low-energy photon spectrometers (LEPS), which are 
highly sensitive to low-energy photons down to a few keV due to a very thin beryllium 
entrance window. They were positioned face-to-face with a distance of only 10 mm to each 
other. The targets were mounted directly between the two detectors to obtain a high detec- 
tion efficiency Figure [6] shows a spectrum of an irradiated naturally composed neodymium 
target. Due to the excellent energy resolution of the detectors (AE™ BM » 0.5 keV at 
E 7 = 100 keV) the observed 7-transitions could be clearly assigned to the corresponding 
/3-decays. However, due to the high detection efficiency, summing effects of coincident 7-rays 
stemming from the same decay cascade cannot be neglected and have to be carefully cor- 
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FIG. 5: The ratio of the normalization factors A 7) o found for the photon flux intensities at different 
energies E mSLX using the reactions 187 Re(7,re) and 11 B(7, 7'), respectively. Both reactions are in 
good agreement with regard to the experimental error bars and the uncertainty of the photoneutron 
cross section of 187 Re (grey error band) [101 ] . The mean ratio of both approaches is indicated by 
the dashed line. At £? max = 7450 keV and -Emax = 7600 keV the 187 Re(7, n) reaction yield was too 
low to obtain data points behind the collimator. 

rected in the analysis of the reaction yield. For the isotopes studied in this work a maximum 
correction of about 10% was found for the summing of 7-rays and x-rays in the case of the 
electron capture of 153 Gd. 

At the second setup a HPGe detector with 30% efficiency relative to a 3" x 3" Nal de- 
tector was used. Here the targets were mounted at a larger distance of 81 mm in front of 
the detector. Thus, almost no summing effects occurred in the spectra. However, due to 
the much smaller detection efficiency this second setup was only used to verify the results 
obtained at the LEPS setup for a few selected targets. 
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FIG. 6: Typical 7-decay spectrum of a neodymium target irradiated at E ma _ x = 8600 keV for a 
duration of about 13 hours. The spectrum was accumulated over a period of 3 hours at the LEPS 
setup. The strongest 7-transitions following the /3-decays of 147 Nd into 147 Pm and 149 Nd into 
149 Pm, respectively, are indicated. 

To determine the absolute efficiency of our setup, the efficiency of the detectors was sim- 
ulated in detail using Geant4. In addition, the efficiencies at certain energies between 
14 keV and 1350 keV were measured using standard calibration sources, which served as 
normalization for the calculated efficiencies. Moreover, a non-calibrated 190 Ir source was 
used to confirm the predicted energy dependence of the efficiency. Figure [7] shows that the 
normalized simulations are in good agreement with the experimental data points. From the 
uncertainties of the experimental data points and the deviation of the normalized simulation 
to these data points, it was estimated that the energy- dependent detection efficiency can be 
determined with an uncertainty of better than 7% for energies of up to about 1500 keV. 
The simulations were also used to account for the self-absorption of the low-energy decay 
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7-rays within the target. Whereas the corrections were almost negligible for the very thin 
neodymium and gadolinium foils, corrections of up to 80% had to be applied for the rel- 
atively thick Sm 2 3 targets. To test the reliability of the simulations the self-absorption 
was also measured at different photon energies for a variety of targets. Deviations between 
measurement and simulation were found to be of the order of a few percent. 
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FIG. 7: Absolute detection efficiency of the HPGe and LEPS setup. Standard calibration sources 
were used to normalize the efficiency curve calculated with Geant4. Due to large summing effects 
of coincident 7-rays at small distances between detector and source, only those calibration sources 
without 7-cascades were used for the normalization in case of the LEPS setup. 
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III. DATA ANALYSIS 



A. Cross section determination 

The reaction yield Y^™* is directly proportional to the energy-integrated cross section 
I a , which is given by the integral from Eq. ([1]). Since bremsstrahlung is characterized 
by a continuous spectral distribution, a deconvolution of the determined integrated cross 
section l a is not possible in general. Therefore, the cross section a( 7)fl ) cannot be derived 
directly. However, if a theoretical prediction for cx( 7jn ) is adopted, one can calculate I a from 
Eq. ([I]) and then derive a normalization factor / for the prediction from a comparison to 
the experimentally determined yields: 



jExp 

/(-Emax) = r B max , : Theory/ P x 7^ ' ( 4 ) 

Js„ ^7^, ^max) • (T^ \E) dE 

By irradiating at different energies E max and deriving / for each energy one can test 
the theoretical prediction for cr( 7 , n ) within different energy ranges. If 0"( 7 , n ) is accurately 
described, then / should be independent of E max an d close to unity. However, our pho- 
toactivation experiments using bremsstrahlung only have limited sensitivity to the shape of 
<7( 7jn ), since the normalization is an average over a wide energy range as illustrated in Fig. 



B. Ground-state reaction rates 



The photodisintegration reaction rate A for a nucleus in a thermal photon bath is given 

by 

POO 

A(T)=/ n^ nck (E,T)a {l>n) (E)dE, (5) 

where 0"( 7>n ) is the photoneutron cross section and ri 7 lanck the photon flux per energy 
interval given by the Planck distribution 

Although it is not possible to produce a thermal photon bath at p-process conditions, 
i. e. with the intensities resulting from temperatures between 2 and 3 x 10 9 K, and thus 
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FIG. 8: The top panel shows the simulated spectral distribution of bremsstrahlung at different 
energies £ max of (1) 8050 keV, (2) 8600 keV, (3) 9050 keV and (4) 9650 keV (solid lines) in 
comparison to the photoneutron cross section cr( 7 n ) of 154 Sm (dashed line) predicted by the NON- 
SMOKER^ 5 code [20(]. The product of ov~ n ) and n 7 yields the integrand of Eq. (JT]) , which covers 



a broad energy range close to the neutron separation energy as shown in the bottom panel. 

to measure stellar reaction rates in the laboratory, we can use two different approaches to 
determine the ground-state reaction rates A g s ' in the experiment. The obvious way is to use 
the normalized cross section cr( 7 , n ) = / • &7 7 ^ ry from the analysis described in the preceding 
section and then calculate the integral of Eq. (j5J). However, using this method one has to 
rely on the adopted theoretical prediction of the energy dependence of the cross section. 
This might lead to large systematic uncertainties, if the predicted shape of the cross section 
deviates significantly from its real shape, e. g. if some resonances above the neutron emission 
threshold are omitted in the theoretical description. Hence, an approach is preferred where 
the reaction rates can be directly determined from the experimental data without the need 
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of any theoretical input. This can be achieved by approximating the Planck spectrum at 
temperature T with a superposition of several bremsstrahlung spectra vS Teras (E, at 



different energies E max [21 ] 



Planck 



(E,T)^^a,(T)n^ ms (E,^ ax ), (7) 



where ai(T) are temperature-dependent weighting coefficients. With this approximation 
Eq. (jSJ) can then be written as 

A(T) a J>(T) /* ^ cms (E,E niax )cT (7 , n) ( J E)^ (8) 

i J 

= I>eoC- ( 9 ) 

i 

Since the integrated cross sections /f* p are directly determined from the experiment the 
reaction rates can be obtained without further assumptions on the energy dependence of 
the cross section. Therefore, this analysis is free of systematic uncertainties stemming from 
the uncertainties of any cross section prediction. The deviation between the approximated 
and the real Planck spectrum in the relevant energy region for astrophysical studies close 
to the neutron threshold energy (the so-called Gamow-like window [211 ]) is of the order of 
10%, depending on how many bremsstrahlung spectra are used for the approximation (see 

Fig. ED. 

IV. RESULTS 

A. Normalization of theoretical predictions 

We have determined normalization factors for two different theoretical predictions of the 
photoneutron cross section cT( 7jn ), which were calculated with the NON-SMOKER WEB code 



by Rauscher 20] and the TALYS code by Koning et al. [22j. Both theoretical predictions 
are based on the Hauser-Feshbach formalism presuming that the nuclear level density in 
the energy range close to the neutron threshold is high enough for a statistical treatment. 
Different results between the two calculations can mainly be ascribed to the nuclear physics 
input used in the codes, e. g. the neutron optical potential, nuclear level densities and the 
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FIG. 9: The top panel shows a thermal Planck spectrum at T = 2.5 x 10 9 K. The weighted sum of 
different bremsstrahlung spectra yields a good approximation of the Planck spectrum within the 
Gamow-like window. The bottom panel shows that this e nerg y region is defined by the product of 



photon flux rKy and photoneutron cross section o"( 7 n ) [211 ] . In this figure cr{-y, n ) was calculated 
with the NON-SMOKER WEB code U for the case of 154 Sm. 



7-ray strength function. 

To compare the results of both predictions we used the default input parameters recom- 
mended by the authors of the codes. These input parameters are mainly intended to provide 
a global description for a wide range of isotopes. The NON-SMOKER WEB code involves the 
neutron optical potential published by Jeukenne et al. 23] with a low-energy modification 
by Lejeune 2^]. The 7-ray strength function is based on a description of Thielemann and 
Arnould 25| using experimental GDR energies and widths if available and the low-energy 
modification of the GDR Lorentzian by McCullagh et al. [26|]. For the nuclear level density 



a global parametrization within the back-shifted Fermi-gas formalism by Rauscher et al. 
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is applied. The TALYS code uses the neutron optical model potential parameterizations of 



Koning and Delaroche 



28j . The 7-ray strength function is obtained from the compilation by 



Kopecky and Uhl 29( and the nuclear level density is also based on an approach using the 
Fermi-gas model |30j. Both the TALYS code and the NON-SMOKER WEB code employ the 
Constant Temperature Model from Gilbert and Cameron 31] for the nuclear level density 
to avoid the divergence of the Fermi-gas model at low excitation energies. 
The experimentally determined normalization factors / of these calculations for different 
energies E max are presented in Fig. [10] and the results are summarized in Tab. |TTJ Both 
statistical model codes are in fair agreement with our experimental data points. A x 2 -test 
confirmed that the normalization factors derived for each isotope can be assumed to be 
independent of the energy E max within the experimental uncertainties. For the absolute 
values of the cross sections we found a mean deviation of 24% and 27% between theory and 
experiment and a mean normalization factor / of 0.82 and 0.73 for the NON-SMOKER WEB 
and the TALYS calculation, respectively. Thus, the predictions seem to slightly overestimate 
the photoneutron cross sections in the rare earth region. 

The quoted uncertainties A/ Sys and A /fit denote the systematic and statistical uncertainties 
of the measurement, respectively. They are discussed in more detail in Sec. IIV CI A/g t was 
derived from averaging the normalization factors / measured at different energies E ma-X with 
regard to the statistical uncertainties of each individual data point. The total uncertainty 
A /total is calculated by a Gaussian error propagation of A/ Sys and A/ fit . Although the ex- 
perimental uncertainty is about 20%, it has to be pointed out that most of the isotopes were 
simultaneously irradiated. Hence, systematic uncertainties are significantly reduced when 
comparing the normalization factors of these isotopes relative to each other. Therefore, the 
discrepancy between experimental data and theory cannot only be due to experimental un- 
certainties. 

Figure [11] shows results from various experiments for the photoneutron cross section in 
the energy region of the giant dipole resonance. These data have been obtained by us- 
ing several experimental techniques like absorption and activation measurements using 
bremsstrahlung, but also direct measurements of the energy dependence of the cross sec- 
tion using a quasi-monoenergetic photon beam produced by the annihilation in fligh t of 



monoenergetic positrons. Details of these experiments are given in 
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33 



34 



35 



36, 3 



The theoretical predictions for the photoneutron cross sections normalized with the factors 
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found in our experiment have been compared to these data. For the isotopes 148 Nd, 154 Sm, 

and 160 Gd the normalized calculations appear to be slightly below the experimental data 

points, but are still fully consistent with these data points within the quoted uncertainties of 

the derived normalization factors. Furthermore, larger deviations have been found for 150 Nd 

and 154 Gd. Unfortunately, in the case of 150 Nd a comparison close to the neutron emission 

threshold is not possible, since no experimental data is available in this energy region. For 

154 Gd a non-negligible photoneutron cross section was even stated below the neutron sep- 

I I 

aration energy of S n = 8894.8 keV in [34] . This indicates systematic uncertainties of these 
experimental data and might explain the discrepancy to the normalized predictions of this 
work. 



TALYS NON-SMOKER WEB 



Isotope 


/ 


A /total 


A/ fi t 


A/sys 


/ 


A /total 


A/ fi t 


A/ sys 


148 Nd 


0.86 


0.176 


0.052 


0.169 


0.86 


0.176 


0.052 


0.169 


150 Nd 


0.55 


0.111 


0.037 


0.105 


0.62 


0.127 


0.047 


0.118 


154 Gd 


0.55 


0.102 


0.052 


0.088 


0.56 


0.100 


0.045 


0.090 


160 Gd 


0.93 


0.151 


0.039 


0.145 


1.15 


0.186 


0.051 


0.179 


154 Sm 


0.74 


0.119 


0.029 


0.115 


0.89 


0.143 


0.034 


0.139 



TABLE II: Normalization factors derived from the experimental data for the theoretical predictions 
of the photoneutron cross section using the TALYS and NON-SMOKER WEB code. The uncertain- 
ties A/fi t and A/ Sys are discussed in detail in the text. A/ tota j denotes the total experimental 
uncertainty derived from a Gaussian error propagation of A/gt and A/ Sys . 



B. Determination of ground-state reaction rates 

We have derived the (7,n) ground-state reaction rates for 148 > 150 Nd, 154 Gd, and 154 Sm 
at temperatures between 2 and 3 x 10 9 K using the approximation of a thermal Planck 
spectrum as described in Sec. IIIIBI The results are presented in Tab. IIII1 
The experimental uncertainties AA| xp Yie id an d AA|'^ p Approx stem from the experimental un- 
certainty of the reaction yield determination and from the uncertainty of the approximated 
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E m [keV] 
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= 0.68 



Nd 



f = 0.55 

XL = - 84 



"Nd 



f = 0.55 
X 2 , = 1-50 




7600 8000 8400 8800 9200 9600 
E [keV] 



NON-SMOKER*™ 

E n [keV] 
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"Gd 



f = 1.15 
1 , = 0-46 



"Gd 



^0.89 
= 0.60 
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0.5 



0.0 



1.5 



7600 8000 8400 8800 9200 9600 
E [keV] 
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1.0 



0.5 



0.0 



FIG. 10: Normalization factors / derived for two different theoretical predictions of the pho- 
toneutron cross section at various energies E mSLX . The grey error band denotes the experimental 
uncertainty A/ to tai discussed in the text. The triangular data points at E mS/X = 9050 keV for 154 Sm 
are independent remeasurements to test the reproducibility of our experimental results. 
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148 Nd 


150 Nd 


154 Gd 


154 Sm 


T = 2.0 x 10 9 K 










Upper energy limit 


8380 


8440 


9850 


8930 



A Exp 

^ A Exp,Yield 

A A g ' s - 

/x Exp,Approx 

\gs. 
TALYS 



\gs. 
A N.S. 



5.97 • 10~ 3 
1.10 • 10~ 3 
1.64 • 10~ 3 
1.09 • 10~ 2 
1.20 • 10" 2 



3.52 • 10~ 3 
0.64 • 10~ 3 
0.85 • 10~ 3 
9.10 • 10~ 3 
9.44 • 10~ 3 



2.74 • 10~ 6 
0.50 • 10~ 6 
0.48 • 10~ 6 
5.10 • 10~ 6 
5.01 • 10~ 6 



T = 2.5 x 10 9 K 



T = 3.0 x 10 9 K 



4.78 • 10~ 4 
0.85 • 10~ 4 
0.54 • 10~ 4 
7.87 • 10~ 4 
6.53 • 10~ 4 



Upper energy limit 8655 


8715 


10120 




9205 


\gs. 
A Exp 


6.45- 


10 1 


4.09' 


10 1 


1.12 • 


io- 


-l 


7.45 


10° 


A \S- S - 

^ A Exp, Yield 


1.17 • 


10 1 


0.74. 


10 1 


0.21 • 


10~ 


-l 


1.32 


10° 


A A gs - 

/v Exp,Approx 


1.39- 


10 1 


0.99 • 


10 1 


0.34 • 


io- 


-l 


0.46 


10° 


\ g - s - 

TALYS 


8.19- 


10 1 


7.74- 


10 1 


2.07' 


io- 


-l 


1.11 


10 1 


\g-s. 
A N.S. 


8.62 • 


10 1 


7.16 • 


10 1 


2.06 • 


10- 


-l 


9.21 


10° 



Upper energy limit 8945 


9010 


10400 


9505 


A Exp 


2.92 • 


10 4 


1.78 • 


10 4 


1.49 ■ 


10 2 


5.04 


10 3 


^ A Exp,Yicld 


0.54- 


10 4 


0.33 • 


10 4 


0.29 ■ 


10 2 


0.88 ■ 


■ 10 3 


A A gs - 

/v Exp,Approx 


0.99- 


10 4 


0.67- 


10 4 


0.52 • 


10 2 


0.28 


10 3 


\g-s- 
TALYS 


3.31 • 


10 4 


3.25 • 


10 4 


2.58 ■ 


10 2 


6.79 


■ 10 3 


A N.S. 


3.39- 


10 4 


2.91 • 


10 4 


2.58- 


10 2 


5.65 


10 3 



TABLE III: The ground-state reaction rates (in s 1 ) determined experimentally by the approxi- 
mation approach of a thermal Planck spectrum for three different temperatures are in good agree- 
ment with two theoretical predictions using the TALYS and NON-SMOKER WEB code. The energy 
range between the neutron emission threshold energy and the quoted upper energy limit (in keV) 
is estimated to contribute about 99% to the^total ground-state reaction rate. AA|'^ pY i c id an d 
^ A Exp Approx accoun t for the uncertainty of the experimental reaction yield and for the uncer- 
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FIG. 11: Comparison of experimental data for the 



todisintegration experiments 
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34 



35|, 



36 



photoneutron cross section from various pho- 



371 ] and theoretical predictions using the TALYS 



(dashed line) and NON-SMOKER WEB code (solid line). The theoretical calculations were normal- 
ized with the factors presented in Tab. [nj Note that the error bands of the theoretical predictions 
due to the experimental uncertainties of the applied normalization factors have been omitted in 
the graphs. 



thermal Planck spectrum within the Gamow-like window, respectively. As will be discussed 
in Sec. II VC} the experimental data points close to the neutron emission threshold suffer 
from large uncertainties. Hence, we discarded these data points for the analysis to increase 
the reliability of the experimental determination of the ground-state reaction rates, although 
the approximation of the Planck spectrum, therefore, became less accurate. Moreover, for 
the isotopes 148 ' 150 Nd and 154 Gd experimental data points were missing at energies further 
above the neutron separation energy. This also reduced the accuracy of the approximation 
of the Planck spectra at temperatures close to 3 x 10 9 K. Thus, except for the case of 154 Sm 
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the approximation in this experiment was less accurate than stated in Sec. IIIIBI For 160 Gd 
too few data points were available in the energy region of interest to derive ground-state 
reaction rates using the approximation approach. 

For comparison the ground-state reaction rates were calculated from Eq. (jSJ) using the unnor- 
malized theoretical predictions for the photoneutron cross sections discussed in the preceding 
section. As shown in Tab. IHII the calculations are in good agreement with the experimental 
results within the experimental uncertainties and the uncertainties of the approximation. 
Consistent with the results found for the normalization factors of the photoneutron cross 
section the statistical model calculations tend to slightly overestimate the reaction rates for 
the studied isotopes. 

C. Discussion of systematic and statistical uncertainties 

An overview of the various experimental uncertainties is shown in Tab. II VI We distinguish 
between systematic uncertainties A sys and statistical uncertainties A sta t- The latter are 
not correlated between the different experimental runs at various energies E max and can, 
therefore, be reduced by a large number of individual measurements. The contribution of 
the various uncertainties to the overall uncertainty will be discussed in detail in the following. 

The dominant uncertainty in our experiment comes from the determination of the photon 
flux as discussed in Sec. Ill Al The systematic uncertainty A sys denotes the mean deviation 
between the normalized Geant4 simulation and the experimental data points stemming 
from the photoscattering reaction 11 B(7,7 / ). Therefore, A sys describes the uncertainty of 
the determined spectral distribution of the photon flux. In addition A sta t accounts for the 
statistical uncertainties of the 11 B(7,7 / ) reaction yields used for the normalization of the 
simulation. 

Close to the neutron emission threshold the reaction yield is highly dependent on the 
maximum photon energy -E max of the activation. Thus, small uncertainties of E mSuX close 
to the neutron emission threshold give rise to large uncertainties in the cross section 
determination as illustrated in Fig. [121 For our analysis we estimated E max to be known 
with an uncertainty of 25 keV. We assumed that there was no systematic deviation of E m&x 
inherent to all experimental runs and, hence, no systematic uncertainty was taken into 
account. Therefore, the uncertainty of E max was only treated statistically in the analysis. 
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Asys 


A s tat 


photon flux 


16% 


5% 


E 


- 


* 


detection efficiency 


4% 


2%t / 15%* 


self-absorption 


5% 




counting 


- 


2% 


7-intensity 


3% 




target mass 




2% 


total 


18% 


6%t / 16%* 



TABLE IV: Systematic and statistical uncertainties of the determination of the experimental (7, n) 
reaction yield. The two uncertainties of the detection efficiency stated refer to the HPGe setup 
(') and the LEPS setup (*), respectively. The uncertainties stemming from an uncertainty in the 
maximum photon energy E ma _ x (*) are illustrated in Fig. [12j The total uncertainty is calculated 
by a Gaussian error propagation. 

Due to averaging over many individual measurements at various E m3uX , the uncertainty 
of E max only represents a minor contribution of the order of a few percent to the total 
uncertainty of the determined normalization factors. 

To derive the reaction yield from the measured activities of the produced unstable isotopes 
the detection efficiency has to be known accurately. As already discussed in Sec. Ill Bl the 
detection efficiency was determined with an uncertainty of 7%. Since each reaction yield was 
derived from several 7-transitions of different energies the average systematic uncertainty 
was reduced to about 4%. Besides, it was found that the observed activation count rate very 
sensitively depended on a proper placement of the target in front of the detectors. Repeated 
measurements with calibration sources showed that an additional uncertainty of 2% at the 
HPGe setup and, due to the very short distance between target and detector, an uncertainty 
of 15% at the LEPS setup had to be taken into account. These measurements also proved 
that this uncertainty could be treated statistically, and, hence, was reduced by averaging 
over many data points. 
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FIG. 12: The expected (7, n) reaction yield of a 154 Sm target (top panel) using a NON- 
SMOKER WEB calculation for the photoneutron cross section (solid line). The dashed and dotted 
lines indicate the error bands assuming uncertainties of E max of 25 keV and 50 keV, respectively. 
The bottom panel shows that close to the neutron threshold this uncertainty clearly dominates the 
uncertainty stemming from the photon flux determination (grey band). 

Other uncertainties stemming from the self-absorption of emitted photons within the target, 
the statistical uncertainty of the counting, the 7-intensities of the decays and the target mass 
were typical of the order of a few percent. 

The total experimental uncertainty derived from the various contributions discussed in this 
section is stated in Tab. IIVI As can be seen, the accuracy is mainly limited by systematic 
uncertainties. Nevertheless, one also has to account for large statistical uncertainties, when 
the reaction yield is determined for a single experimental run. However, the statistical 
uncertainties are significantly reduced when averaging over many data points in the analysis 
as done for the deduction of the normalization factors and ground-state reactions rates. 
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Therefore, the total experimental uncertainties for the results presented in this work are 
typical of the order of 20%. 

V. SUMMARY AND CONCLUSION 

In the last years extensive studies of photoneutron reactions with astrophysical implica- 
tion for isotopes in the mass region A > 186 have been carried out, but no experimental data 
have been provided for the rare earth isotopes so far. However, reliable experimental data 
for a wide range of isotopes in various mass regions are mandatory in order to constrain and 
to improve the theoretical calculations required for astrophysical reaction networks. The in- 
tention of this work was to test the theoretical predictions of two different statistical model 
codes for the photoneutron cross sections and ground-state reaction rates in the rare earth 
region. Although a variety of photoneutron reactions in this mass region was studied in the 
70s, these experiments mainly focused on the investigation of the Giant Dipole Resonance. 
Hence, data points stemming from these experiments usually underlie large statistical and 
systematical uncertainties close to the neutron threshold and do not serve as a reliable test 
for theoretical predictions in the astrophysically relevant energy region. Therefore, we de- 
termined normalization factors for the calculations of the photoneutron cross sections of the 
isotopes 148 ' 150 Nd, 154 Sm, and 154 ' 160 Gd in photoactivation experiments close to the neutron 
separation energy. Moreover, we derived (7, n) ground-state reaction rates in the astro- 
physically relevant energy region for the p process using a superposition of bremsstrahlung 
spectra of various energies. 

The need of experimental data of photoneutron reactions for astrophysical reaction networks 
is twofold. First, for ^-process studies reliable predictions of reaction rates for a wide range 
of isotopes need to be provided. Since most of the isotopes involved in the j9-process reac- 
tion network are unstable, the nuclear physics input for these isotopes is rare and, hence, 
theoretical calculations usually cannot be adjusted locally. Therefore, the reliability of the 
theoretical models should be tested using only global input parameters. In this work we 
found that in the rare earth region the predicted cross sections and the reaction rates for 
the selected isotopes agree with the experimental data within a factor of two. On average, a 
mean uncertainty of less than 30% has been determined. This is consistent with the results 
found in previous experiments in the mass region A > 186. In addition, no systematic 
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dependence of the deviation with atomic mass has been observed. On the one hand, this 
proves the reliability of the statistical model codes for a wide range of isotopes along the 
valley of stability. But it also indicates an appropriate predictive power for the reaction 
rates of unstable isotopes, since only global input parameters were used in the calculations . 
However, for almost all studied isotopes the cross sections were found to be overestimated by 
theory, which calls for further investigations. Moreover, it is planned to study the isotopic 
chain of the cerium isotopes in the near future to test the validity of the statistical model 
in this mass region for isotopes close to and at the neutron shell closure. Although the 
predictions agree fairly well with the experimental data of this work, it needs to be empha- 
sized that experimental studies cannot account for the stellar enhancement of cross sections 
and reaction rates due to the thermal population of low-lying levels in the parent nuclei 
under stellar conditions. This so-called stellar enhancement factor can only be derived from 
theory. 

The second astrophysical application of experimental studies of photoneutron reactions arises 
in the investigation of branching points within the s process. In this context the photoneu- 
tron reactions of 148 Nd and 154 Gd are of particular interest among the studied isotopes of 
this work. Since the neighbouring branching points 147 Nd and 153 Gd cannot be studied in 
neutron capture experiments due to their short half-life, the aim is to derive the neutron cap- 
ture cross sections of these branching points from the photoneutron reaction of the isotopes 
148 Nd and 154 Gd, respectively. It was shown in this work, that the two adopted theoretical 
predictions overestimate the photoneutron cross sections by up to a factor of two. This 
might indicate that also the predicted neutron capture cross section of the inverse reactions 
need to be adjusted accordingly when being calculated within the same statistical model 
code. For astrophysical studies, a direct correlation between the stellar neutron capture and 



photodisintegration rate is described by the so-called principle of detailed balance [38|, 139 ]. 
This correlation only holds under the assumption that the low-lying levels in both the parent 
and the residual nucleus are thermally populated to a significant amount. This condition is 
fulfilled under sufficiently hot temperatures in a stellar environment. However, experiments 
in the laboratory only address transitions starting from the ground state of a certain nu- 
cleus and, therefore, the principle of detailed balance is not applicable any more for a direct 
comparison of the laboratory neutron capture and photodisintegration rate. Hence, without 
any further assumptions, conclusions for the neutron capture rate can hardly be drawn from 
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the reaction rates derived in photodisintegration experiments. 

Instead of directly deriving the neutron capture reaction rate from experiment a promising 
approach is to improve the nuclear physics input of the statistical model codes in order 
to increase the reliability of the theoretical predictions. To obtain the most accurate cal- 
culations for a single nucleus the input parameters should be adjusted locally. Therefore, 
the aim of future studies should be to provide improved nuclear physics input for 148 Nd 
and 154 Gd, while the experimental results of this work serve as a reliable test for any new 
complete set of input parameters. It can be assumed that a set of input parameters yielding 
improved predictions for the photoneutron reaction, will then also enhance the calculations 
of the neutron capture cross sections of the branching points 147 Nd and 153 Gd. At this 
point, however, we would like to make clear that the experimental data provided by this 
work set constraints on the absolute value of the photoneutron cross section, but do not 
allow for separately adjusting the various nuclear physics parameters such as the neutron 
optical potential, nuclear level densities and the 7-ray strength function. Therefore, further 
experimental investigations with direct access to these parameters are mandatory in the 
future. 
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